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ABSTRACT 


This  report  presents  additional  theoretical  results  based  on  a  visco¬ 
plastic  flow  model  for  hypervelocity  perforation  of  thin  plates^  which  was 
given  in  References  1  and  2.  The  solutions  of  the  governing  equations  in¬ 
cluding  the  yield  strength  are  obtained.  The  results  of  a  parametric  study 
performed  on  an  IBM  1620  computer  are  also  included.  The  effects  on  the.rai. 
dius  of  perforation  due  to  various  physical  quantities  are  discussed.  The 
theoretical  radius  of  perforation  and  residual  velocity  are  compared- with 
the  experimental  values  -  obtained  by  R,  W.  Watson. 


I.  INTRODUCTION 


In  references  1  and  2,  a  visco-plastic  flow  theory  has  been  pre¬ 
sented  for  the  solution  of  the  hypervelocity  perforation  of  thin  plates. 
In  this  approach,  the  material  is  assumed  to  behave  like  an  incompress¬ 
ible  visco-plastic  material.  The  perforation  process  is  described  by  a 
simple  shear  plugging  model.  The  only  non-vanishing  strain  component 
is  the  shear  strain,  which  is  assumed  constant  across  the  thickness,  of 
the  plate.  The  resulting  equation  of  motion  is  of  the  diffusion  type 
and  is  solved' by  the  Laplace  transform  methods.  The  diameter  of  the 
perforated  hole  in  the  plate  is  then  predicted  by  a  separation  criteri¬ 
on.  For  simplicity,  the  yield  stress  k  was  dropped  from  the  equation. 

In  the  present  report,  the  solutions  of  the  basic  equations,  in¬ 
cluding  the  yield  strength,  have  been  obtained  and  expressed  in  infin¬ 
ite  series  forms.  A  parametric  study  is  then  presented,  in  which  sev¬ 
eral  nuTiierical  values  are  assigned  to  each  physical  parameter  and  their 
effects  on  the  radius  of  perforation  determined. 


II.  BASIC  EQUATIONS 


Under  the  assumption  of  the  simplified  perforation  model,  the  only 
equation  of  motion  is 

(1) 

The  visco-plastic  property  of  the  material  is  assumed  to  be 

The  relations  between  the  strain  ,  axial  displacement  and  the  velocity 
arc 
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Combining  Equations  1,  2  and  3  yields  the  following  governing  dif¬ 
ferential  equation, 


r  ar  r  p. 

The  initial  and  boundary  conditions  raay  be  expressed  as 


u=  o 

AT 

t ' r>a 

U  -  Ur 

AT 

ri¬ 

ll 

^i 

-J 

tl 

H.>jl  “Kti- 

II 

O 

t>  O  5  1^  =  61 

u— ^  o 

AT 

Y'— ®o 

wheye 

k’. 

Equations  4  and  5  may  be  made  dimensionless  by  introducing  the  Rey¬ 
nold's  number  and  Bingham's  number, 

R  5  Reynold's  Number  =  — — 


B  r  Bingham's  Number 


and  the  dimensionless  quantities:  f’  ,  and  v  defined  by 

r'=  ;  t  =  ;  u  = 


a  ’ 

The  resulting  equations  are 


t'  =  t 
^  a. 


U"'  s  O 
ir'  =  I 

-§^-8-Rf^  =  0 

1/'  =  o 


AT  t  >Oj  r  »  ) 
AT  t*  =  Oj  r'>  I 
AT  t'  =  o,  r  '=  j 
AT  t'?-o^r'*i 
AT  tVo^r’««» 


where  H=  ^  dimensionless  mass  ratio. 


Equations  6  to  10,  with  BrO,  have  previously  been  solved  and  the 
results  reported  in  Reference  2.  In  this  report,  solutions  are  obtained 
for  the  case  B/O  by  using  Laplace  transform  and  asyraptotic  expansion 
techniques.  The  results,  which  are  summarized  in  Appendix  I,  are  in  the 
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forra  of  equations  for  tr*.,  Z.  ,  in  terms  of  t  and  T 

For  any  given  set  of  values  of  parameters  R,  B  and  H,  these  dimension¬ 
less  quantities  may  be  calculated  from  Equations  1-14,  1-15,  1-16  and 
1-17,  respectively.  The  calculations  involved  have  been  programmed  on 
an  IBM  1620  computer.  The  output  from  the  computer  is  either  in  tab- 

/  t 

ular  form  or  in  graph  form,  plotted  as  vs.  r'  and  -^—7 vs.  r' 

curves  at  different  time  intervals.  A  typical  set  of  curves,  plottfed 
by  the  computer,  are  shovm  in  Figures  1  and  2. 

The  same  separation  criterion  as  presented  in  Reference  2  is  used 
here  for  the  prediction  of  the  perforation  diameter.  This  criterion 
stipulates  that  separation  occurs  when  the  plate  material  changes  from 
a  "viscous  fluid"  to  a  "non-viscous"  fluid.  It  is  further  assumed  that 
the  material  is  viscous  when  the  viscous  stress  is  larger  than  the  yield 
stress;  the  material  is  considered  a  fluid  only  if  the  strain  is  larger 
than  a  constant,  in  this  case  the  static  ultimate  shear  strain  is  used. 

In  equation  form,  this  separation  criterion  states  that  separation  occurs 
when 

^  ^  (11) 

The  actual  procedure  in  applying  this  criterion  and  in  finding  the 
radius  of  perforation  is  sho’im  in  Figures  1,  2  and  3.  In  Figure  1,  a 
horizontal  dotted  line  is  drawn  for  a  particular  value  of  From 

the  intersections  between  this  line  and  the  strain  vs.  radius  curves, 
the  critical  strain  curve  may  be  constructed  as  shovm  in  Figure  3.  Sim¬ 
ilarly,  the  critical  strain-rate  curve  may  be  constructed  from  Figure  2. 
The  intersection  between  the  critical  strain  and  critical  strain-rate 
curves  gives  the  radius  of  perforation  Tp  and  the  timetp  when  separation 


occurs . 
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III,  PAi’vAl-IETaiC  STUDY 

A  parametric  study  was  performed  in  order  to  find  out  the  effect 
on  the  perforation  process  due  to  different  parauiecers .  The  parameters 
used  were  dimensional  physical  quantifies,  such  as  plate  thickness  (inch.) 
and  projectile  velocity  (ft/sec),  rather  than  dimensionless  quantities 
B,  R,  and  H.  This  approach  was  used  to  obtain  directly  the  effect  due 
to  physical  quantities.  The  follov7ing  parameters  V7ere  studied:  projec¬ 
tile  velocity  Vo>  projectile  diameter  a,  projectile,  mass  plate  thick¬ 
ness  h,  plate  density  9  ,  yield  strength  k,  coefficient  of  viscosity ja  , 
ultimate  shear  strength  ^  and  ultimate  shear  strain  A)  .  A  set  of  base 
values  for  these  parameters  was  first  chosen,  the  strain  and  strain-rate 
distributions  were  calculated  and  the  radius  of  perforation  determined. 
Then,  each  parameter  was  varied  individually  with  all  the  other  para¬ 
meters  held  at  their  base  values.  The  effect  on  the  radius  of  perfor¬ 
ation  due  to  each  parameter  was  evaluated,  and  discussed. 

The  set  of  base  values  used  are 

Vo  =  3170.  ^/sec  ji.  ■  100  Ib-sec/ft^ 

a  ■  1/32"  k  r  24,000  psi 

h  r  1/16",  1/8",  3/16",  1/4"  9  =  5.3707963  slug/ft^ 

Ml  r  .025  gm  TT,  ;  91,000  psi 

(To  =  .02 

The  values  for  Vq,  a,  h  and  Mj^,  are  similar  to  those  used  in  a  series 
of  experiments  performed  by  Watson,  as  shown  in  Reference  3.  The  values 
for  the  rest  of  the  parameters,  all  pertaining  to  nuterial  properties, 
were  either  taken  from  handbooks  or  estimated.  A  coefficient  of  viscos¬ 
ity  of  100  lbs-sec/ft2  ig  believed  to  be  a.  reasonable  estimated  value  as 
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discussed  in  Reference  1.  The  yield  strength  in  shear  of  the  plate  and 
the  density  of  the  plate  are  those  properties  of  2024-T3  aluminum.  An 
ultimate  strength  of  the  plate  of  91,000  psi  was  assumed  for  dynamic 
conditions.  Finally,  the  critical  strain  of  .02  was  used  because  the 
ultimate  strain  of  2024-T3  aluminum  under  static  conditions  is  of  the 
order,  of  magnitude  of  one  per  cent. 

IV.  RESULTS  AND  DISCUSSIONS 

The  numerical  values  used  for  the  parameters,  together  with  the 
calculated  results,  are  presented  in  Table  1.  These  results  are  also 
plotted  in  curves  as  shown  in  Figures  4  to  10.  In  these  figures,  the 
radius  of  perforation  r  p  jg  plotted  as  ordinate  and  the  plate  thickness 
h  as  abscissa,  with  one  particular  physical  quantity  as  parameter.  Three 
curves  are  shown  on  each  figure  representing  three  values  of  the  para¬ 
meter  under  consideration,  one  at  the  base  value,  one  above  and  one  be¬ 
low.  These  curves  are  discussed  in  the  following  paragraphs.  In  ad¬ 
dition,  the  combined  effect  of  the  ultimate  strength  and  ultimate  strain 
are  discussed  in  Section  10.  Comparisons  with  test  data  are  given  in 
Sections  11  and  12. 

1,  Plate  Thickness 

In  Figures  4  Co  10,  it  can  be  seen  that  an  increase  in  plate 
thickness  decreases  the  r'p,  regardless  of  the  other  parameters  in¬ 
volved. 

2.  Plate  Density 

'  Figure  4  shows  the  effect  of  the  variation  of  the  plate  density 
on  the  dimensionless  perforation  radius.  By  increasing  the  density 


-6 


from  5.37  slug/ft^  to  10.74  slug/ft^  the  r'p  is  decreased  through¬ 
out  the  four  thicknesses.  A  decrease  in  the  density  to  2.685  slug/ft^ 
increases  r'p  over  all  plate  thicknesses.  This  range  of  density  in¬ 
cludes  such  materials  as  aluminum,  beryllium  and  magnesium, 

3.  Coefficent  of  Viscosity 

The  effect  on  r'p  due  to  the  variations  of  the  coefficent  of 
viscosity  is  shown  in  Figure  5.  A  decrease  in  r!p  is  caused  by  a 
decrease  in  viscosity  from  100  Ib-sec/ft^  to  50  Ib-sec/ft^  and  an 
increase  is  observed  when  the  viscosity  was  raised  to  200  Ib-sec/ft^. 
This  range  was  found  to  be  sufficient  to  evalutate  the  trend.  The 
larger  value  of  viscosity  is  seen  to  have  a  considerable  effect  on 
the  r'p. 

4.  Projectile  Velocity 

Figure  6  shows  that  an  increase  of  the  projectile's  initial 
velocity  from  10,400  fps.  to  20,800  fps  enlarges  the  r'p  while  a 
decrease  to  7800  fps  decreases  the  r'p.  This  range  of  velocities 
represents  the  current  laboratory  capability  aiid  perforation  test 
data  at  this  velocity  range  are  expected  to  be  available  shortly. 

5 .  Ultimate  Strength 

In  Figure  7,  the  ultimate  strength  (  )  is  decreased  from 

91,000  psi  to  65,000  psi  resulting  in  an  increase  of  r'p.  With  Xfi 
increased  to  125,000  psi  the  r'p  is  decreased.  This  range  was  used 
due  to  the  uncertainity  of  the  dynamic  value  of  the  ultimate  strength. 

6 .  Critical  Strain 


The  critical  strain  (  Jfo)  effect  on  r'p  is  shown  in  Figure  8. 
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The  r'p  Is  increased  for  an  increase  of  from  .02  to  .035.  For 
a  decreased  value  of  Zf*  of  .005  the  r'p  is  decreased.  This  range 
of  jC,  represents  a  reasonable  variation  of  . strain  from  the  static 
ultimate  strain  of  .01. 

7.  Projectile  Mass 

In  Figure  9  the  effect  of  the  variation  of  projectile  mass  on 
the  dimensionless  perforation  radius  is  shown.  When  the  mass  is  in- 

"6  —A 

creased  from  1.71  x  10  slugs  to  3.42  x  10  °  slugs  a  larger  r'p  is 
observed.  With  a  decrease  of  the  mass  to  .86  x  10"^  slugs  a  smaller 
r'p  is  obtained.  The  radius  of  the  projectile  was  kept  constant 
for  these  three  values  of  the  projectile  mass.  This  change  of  mass 
may  be  accomplished  by  a  change  of  either  length  or  density  of  the 
projectile. 

8.  Projectile  Radius 

The  variation  of  rp  with  the  change  in  projectile  radius  is 
shoum  in  Figure  10.  Since  r'psrp/a^  it  is  necessary  to  show  the 
dimensional  rp  variation  with  a.  An  increase  in  the  projectile's 
radius  from  .03125  inch  to  .0442  inch  produces  a  slight  increase  in 
tp.  While  a  decrease  to  .0221  inch  causes  a  small  decrease  in  rp. 
However,  the  differences  between  all  three  are  rather  small.  The 
mass  of  the  projectile  was  held  constant-  throughout  this  series, 
therefore,  this  change  in  the  projectile's  radius  would  be  accomp¬ 
anied  by  a  change  in  either  the  length  or  density  of  the  projectile. 

9.  Yield  strength 

When  k  is  decreased  from  24,000  psi  to  5000  psi  the  r'  is  un- 

ir 

changed.  With  the  exception  of  one  case  (Run  #39,  h  s  1/4"),  in¬ 
creasing  k  to  43,000  psi  does  not  change  the  perforation  radius. 
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For  Run  #39,  the  r'p  is  slightly  increased.  This  range  of  yield 
strength  values  was  vsed  to  account  for  the  possible  variation  due 
to  the  dynamic  effect. 

10.  Combined  Effect  of  ifo  and 

In  Figure  11  it  can  be  seen  that  a  particular  diameter  can  be 
obtained  from  the  critical  strain  and  strain-rate  curves  by  varying 
either  the  critical  strain  and/or  ultimate  strength  of  the  plate 
material . 

a)  Critical  Strain  Constant 

An  increase  in  ultimate  strength  from  91,000  psi  to  125,000 
psi  will  decrease  both  the  r'p  and  t'p.  Decreasing  the  ul¬ 
timate  .strength  to  65,000  psi  increased  the  r'p  and  t'p. 

b)  Ultimate  Strength  Constant 

Increasing  critical  strain  from  .02  to  .035  increased  the 

t'p  and  t'p.  A  decrease  to  .005  decreased  these  same  quan¬ 
tities. 

c)  Variation  of  Both  Strain  and  Ultimate  Strength 

By  observing  Figure  11  one  can  see  that  many  combinations 
of  critical  strain  and  ultimate  strength  yield  a  particur 
lar  value  for  t'p.  However,  the  t'p  varied  for  each  com¬ 
bination  of  critical  strain  ^d  ultimate  strength. 

11.  Residual  Velocity 

The  residual  velocity,  (or  the  velocity  of  the  projectile  after 
the  perforation  process  is  completed),  is  calculated  for  a  few  cases. 
These  are  shown  in  Figure  12  in  which  the  velocity  tr,  is  also  presented. 
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The  velocity  xTo  is  calculated  by  a  simple  momentum  exchange  between 
the  projectile  and  the  portion  of  the  plate  directly  in  contact  with 
it.  The  quantity  M  (vQ-Vrea)  represents  the  momentum  transferred 
to  the.  perforated  plate  fragments  (spall). 

Comparing  these  curves  with  the  test  data  of  Reference  3,  an 
apparent  discrepancy  appears.  In  particular,  it  may  be  observed  in 
Figure  12,  that  at  a  plate  thickness  of  .31  cm,  the  test  residual 
velocity  is  much  higher  than  the  theoretical.  There  are  two  poss¬ 
ible'  explanations  for  this  discrepancy.  First,  the  spall  may  flash 
back  in  a  direction  opposite  to  that  of  the  projectile,  and  as  a 
result,  the  forward  residual  momentum  and  velocity  are  increased; 
or,  second,  the  portion  of  the  plate  directly  in  contact  with  the 
projectile  may  leave  the  plate  with  a  velocity  less  than -the  res¬ 
idual  velocity  of  the  projectile.  However,  according  to  Watson's 
experiments,  there  was  no  flash-back  observed  in  the  perforation 
of  aluminum  plates.  Presently,  the  Bureau  of  Mines  is  conducting 
tests  to  determine  the  residual  velocity  of  the  spall  in  order  to 
evaluate  the  feasibility  of  the  second  explanation  regarding  this 
inconsistency. 

12.  Comparison  of  r'p  With  Test  Data 

Observation  of  Figure  13  shows  the  comparison  of  the  perfor¬ 
ation  diameter  of  various  thicknesses  of  aluminum  2024-T3  plate. 

At  a  plate  thickness  of  1/16"  there  is  a  discrepancy  with  the  the¬ 
oretical  value  being  slightly  higher.  However,  for  the  other  plate 
thicknesses  (1/8",  3/16"),  there  is  general  agreement. 
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FOR  BASE  VALUES  OF  ULTIMATE  STRENGTH  (T^>) 
AND  CRITICAL  STRAIN  (3o) 
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FIG.  5  THE  EFFECT  OF  THE  COEFFICIENT  OF 

VISCOSITY (^l)  ON  THE  PERFORATION  RADIUS 
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FIG.  8  THE  EFFECT  OF  CRITICAL  FIG.9  THE  EFFECT  OF  THE  MASS 
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FOR  VARIOUS  VALUES  OF  ULTIMATE  STRENGTI^  (T«) 
AND  CRITICAL  STRAIN  (A) 


VELOCITY 


Mis  .025  gm 
Vo=3l70  */,,c 
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CARNEGIE  TEST  RESIDUAl*  VELOCITY 
^-THEORETICAL  Vo  =  --^7 '  . 

M  total 

©-THEORETICAL  RESIDUAL  VELOCITY 


FIG.  12  RESIDUAL  VELOCITY  VS.  PLATE  THICKNESS 
FOR  202  4“T3  ALUMINUM 
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ORIGINAL  VELOCITY  =3170 

RADIUS  OF  PROJECTILE  = 

PLATE  material  =2024-T3  AL 
MASS  OF  projectile  =  .025  GM. 

k  EXIT 


h  -PLATE  THICKNESS 


(CARNEGIE) 

(CARNEGIE) 

(DREXEL) 


FIG.I3  COMPARISON  OF  PERFORATION  DIAMETERS- 
CARNEGIE'S  TEST  RESULTS  VS.  THEORETICAL 
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VIII.  NOTATIONS 


a  z 

6  > 

C  = 

h  a 
H  a 

k  = 

M  a 

Ml  a 

r  a 


r’  s 


r'D  a 

R  a 

t  s 
t'  : 

V  a 

V'  a 
Vo 
z  : 
z'  z 
«.  z 

f  = 

"O  a 

asa;  - 
ap‘  " 


Radius  of  projectile 
Bingham  No.  s  ■ 

Constants 

Thickness  of  plate 


Mass  Ratio 


Yield  Strength  in  Shear 

Ml  -t-  .-rra'^K  f 

Maas  of  projectile 
Radial  distance 

Dimensionless  radial  distance  a 
Dimensionless  perforation  radial  distance 
Reynold's  No.  a 
Time 

Dimensionless  time  s 

“3“ 

Velocity 

Y 

Dimensionless  Velocity  a 
Initial  velocity 
Displacement 

Dimensionless  displacement 
critical  Strain 
Density  of  plate  --  ■' 

Coefficient  of  Viscosity 
Kinematic  Coefficient  of  Viscosity  r 
Dimensionless  Strain  -  rate 


a  Dimensionless  Critical  Strain  -  rate 


Dimensionless  Strain 
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VIII.  APPENDIX  I:  Solution  of  the  Governing  Equations 


Applying  the  Laplace  transform  to  Equation  6  results  in 


.  j_  Su-  _  J _ ^ 

r'  3r‘  r*  p 


u 


(i-l) 


where  tr  is  the  transform  of  u- j  x  Rp  and  p  is  the  parameter  of  the 
transform.  Introducing 

■r  =  <^r‘  ,  a  *  an»o  C=; 

Equation  ( I- 1)  becomes 

(1-2) 


i  au  c 

«--= — u  = 


The  boundary  conditions,  Equations  9  and  10,  transform  into 

q#--  p--  =  o  (1-3) 


AS 


oo 


XT  =  O 


(1-4) 


The  initial  conditions  are  satisfied  by  the  transforms  (1-2)  and 
(1-3),  The  general  solution  of  the  homogeneous  equation  associated 
with  Equation  (1-2)  is 

ij  =  AKoCr)  + 

where  Iq  and  Kq  are  the  modified  Bessel  Functions  of  the  zero  order, 
and  of  the  first  and  second  kind  respectively;  A  and  Bq  are  constants. 
The  particular  solution  for  Equation  (1-2)  is  obtained  by  the  method  of 
variation  of  parameters.  The  Wronskian  of  the  general  solution  of  the 
homogeneous  equation  is 

W  «  <  KolJ  -  loKo  )  =  ^ 


(1-6) 


1 


where  the  prime  represents  differentiation.  It  can.  be  shown  that  the 
particular  solution  of  Equation  (1-2)  is 

C[l,  jK^dlr  -Kojiocif'] 

The  general  solution  of  Equation  (I-.2)  is  thus 

u-=  AKoCr)  ^  +  c[lo5Kodr  -  Kojlod^^  ]  (1-8) 

At  r  =  c«  ,  the  functions  Kq,  Io_yKo<j^  and  Ko|  all  vanish. 

.Therefore,  to  satisfy  Equation  (1-4),  it  is  only  necessary  to  set  BqcO. 
Substituting  Equation  (1-8)  into  Equation  (1-3)  yields  the  constant  A, 

p  K*a)- J± 

and  the  solution  of  Equation  (1-2)  which  satisfies  Equations  (1-3)  and 
(1-4)  is 

0=  AKoCp)  >c[i«  jKodf^-  KoJ^Iod?-] 

Because  only  a  small  time  interval  is  of  interest,  the  following 
asymptotic  expansions  of  Kq  and  Kj.,  are  used  (Reference  4): 

K.Cf)  =  €  [  '  -  ^  ♦ . } 


K,cr)  =  U^)  e  I"  t-v  ^  - . J 

The  resulting  equation,  after  algebraic  simplification,  is 


(1-12) 


T’  -  >  .  ■  1  1 

8r'  '  ^ 

.JL 

fcl-r' 

B  , 

.  -  - 

R.*  r '  * 

Similarly,  expansions  for  lo  and  can  be  obtained. 
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using  Standard  transform  tables  (Reference  5)  gives 


The  corresponding  displacement,  strain  and  strain-rate  are 
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where , 


Cl  « 


C2»  (h^-4h-6H-^9^5) 


Its  r 


:.  * 


-  l)‘r  *  ('Iw  *;554  )f.  ’ 
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C5  « 


3_ 
8P  ' 


C7= 

Co-  (h*-  tH*-  H  SH  B 

'  ('  i»^‘‘  B  ‘•xj^  )  7T-.»  * 


.  4—  w  ■  *o^  ^~x  - _4lli_ 
Vi024^  /r*^  TvtfS'p’'* 


